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ABSTRACT 

Context. Oxygen sequence Wolf-Rayet (WO) stars represent a very rare stage in the evolution of massive stars. Their 
spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The 
properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive 
stars, as well as their remaining lifetime and the expected properties of their supernovae. 

Aims. We aim to homogeneously analyse the currently known presumed-single WO stars to obtain the key stellar and 
outflow properties and to constrain their evolutionary state. 

Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code CMFGEN to model the X-Shooter 
spectra of the WO stars and deduce the atmospheric parameters. We calculate dedicated evolutionary models to 
determine the evolutionary state of the stars. 

Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface 
helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models 
with helium zero-age main sequence masses of Mne.ini = 15 — 25Mq that exhibit a fairly strong (on the order of a few 
times lO“®M 0 yr“^), homogeneous (/c > 0.3) stellar wind. 

Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of Mini = 40 —OOMq. 
They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years. 


1. Introduction 

The enigmatic oxygen sequence Wolf-Rayet (WO) stars 
represent a very rare stage in massive star evolution. Their 
spectra are characterized by strong emission lines of highly 
ionized carbon and oxygen, and in particular a strong O vi 
A3811-34 A emission line. Their emission-line spectra point 
to dense, outflowing atmospheres. Despite their rarity, the 
WO stars can provide key information in our understanding 
of massive star evolution. 

Although the progenitors of hydrogen-free type Ib and 
Ic supernovae (SNe) have not yet been identified, it is very 
likely that they are evolved Wolf-Rayet stars (e.g., Yoon 
et al. 2012). In particular, WO stars are potential progeni¬ 
tors of the helium-deficient type Ic SNe, as they may have 
a very low helium abundance. If they are rapidly rotating 
when they explode, they could produce an associated long- 
duration gamma-ray burst (e.g., Woosley & Bloom 2006). 

The tell-tale signature of WO stars is their O vi A3811- 
34 A emission. This emission was first found in the spec¬ 
tra of the central stars of planetary nebulae. However, 
Sanduleak (1971) pointed out that five of the stars showing 

* Based on observations obtained at the European Southern 
Observatory under program IDs 091.G-0934 and 093.D-0591. 


Ovi A3811-34 A have broad, Wolf-Rayet like emission lines 
and do not appear to have an associated nebula. It was 
therefore suggested that they are part of the carbon se¬ 
quence Wolf-Rayet (WC) class. Barlow & Hummer (1982) 
argued that these WC O vi stars should be seen as a sepa¬ 
rate class of Wolf-Rayet stars, and introduced the first WO 
classification scheme. 

Since their discovery, the WO class has been commonly 
interpreted as a very short stage of evolution of mas¬ 
sive stars covering an initial mass range of approximately 
45 — 60Mq after the carbon sequence Wolf-Rayet (WC) 
phase (e.g., Sander et al. 2012; Langer 2012). In such a 
scenario, the emission lines of highly ionized oxygen reflect 
the high oxygen abundance that is expected near the end 
of core-helium burning. The very high stellar temperature 
that is needed to produce O vi emission is expected if WOs 
are indeed the descendants of WC stars, as the envelope is 
being stripped by the stellar wind and consecutively hot¬ 
ter layers are revealed. An alternative scenario is that WO 
stars originate from higher mass progenitors. In this case 
the stars are hotter and the O vi emission could purely be 
an excitation effect, and a high oxygen abundance is not 
necessarily implied (e.g., Hillier & Miller 1999). 
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Table 1. Overview of the known WO stars. 


ID 

R.A. 

(J2000) 

Dec. 

(J2000) 

Spectral type" 

Z 

(Zq) 

SIMBAD ID*’ 

Other IDs 

WR102 

17:45:47.56 

-26:10:26.9 

W02 

1 

V* V3893 Sgr 

[S71d] 4 (Sand 4) 

WR142 

20:21:44.35 

-k37:22:30.6 

W02 

1 

WR 142 

[S71d] 5 (Sand 5) 

WR93b 

17:32:03.31 

-35:04:32.4 

W03 

1 

WR 93b 


BAT99-123 

05:39:34.31 

-68:44:09.1 

W03 

0.5 

Brey 93 

[S71d] 2 (Sand 2) 

LH41-1042 

05:18:11.01 

-69:13:11.3 

W04 

0.5 

[L72] LH 41-1042 


LMC195-1'^ 

05:18:10.33 

-69:13:02.5 

W02 

0.5 

— 


DRl 

01:05:01.61 

-k02:04:20.6 

W03 

0.15 

NAME DR 1 in IC 1613 

[BUG2007] B 17 

WR30a 

10:51:38.93 

-60:56:35.2 

W04 -t 05((f)) 

1 

V* V574 Car 

[MS70] 4 (MS4) 

Skl88 

01:21:04.13 

-73:25:03.8 

W04 -1- 04 

0.2 

2MASS J01310412-7325038 

[S71d] 1 (Sand 1) 


Notes. The upper part of the table shows the (apparently) single WO stars, the bottom part the binaries. 

See Section 2.1 for the spectral classihcation of the single stars. Binary classihcations are from Massey et al. (2000) and Moffat 
& Seggewiss (1984). http://simbad.u-strasbg.fr/simbad/ Not yet listed in Simbad. Designation from Massey et al. (2014). 


If WO stars represent a later stage of evolution com¬ 
pared to WC stars, the surface abundance of carbon and 
oxygen is expected to be high. If these can be measured, 
and a good estimate of the stellar mass can be obtained, the 
WO stars can provide unique observational constraints on 
the elusive ^^C(q;, 7)^®0 thermonuclear reaction rate. This 
rate is currently only weakly constrained, with an uncer¬ 
tainty on the order of 30% (e.g., Tur et al. 2006). 

Currently, nine members of the WO class are known, 
two of which are in a binary system. Table 1 lists the co¬ 
ordinates, names, and metallicities (Z) of all these stars. 
Two WO stars have recently been discovered in the LMC. 
The first, the W04 star LH41-I042, was discovered by 
Neugent et al. (2012). The second one was reported by 
Massey et al. (2014), and is, remarkably, located only 9" 
away from LH41-1042. 

As a first step in our effort to investigate the nature 
of the WO stars, Tramper et al. (2013, henceforth Paper 
I) performed a detailed spectroscopic analysis of DRl. 
Located in the Local Group dwarf galaxy IC 1613, this star 
is the lowest metallicity WO star known, with a metallicity 
Z ~ 0.1-0.2 Zo (Paper I). The stellar parameters that 
were derived confirm the expected very high stellar tem¬ 
perature, close to the helium terminal-age main sequence 
(He-TAMS). However, the derived surface abundances of 
helium, oxygen, and carbon are comparable to those found 
in early WC stars. Thus, if DRl is representative for the 
WO stars, these stars may be descendants of higher mass 
progenitors than WC stars. 

In this paper, we perform a homogeneous spectroscopic 
analysis of all the remaining (apparently) single WO stars 
that are known, with the exception of the recently discov¬ 
ered star in the LMC (LMC195-1, Massey et al. 2014). The 
stellar parameters, together with the results from paper 
I, are used to determine the nature of WO stars. We use 
dedicated helium-burning models to constrain their evolu¬ 
tionary stage as well as to predict their remaining lifetime. 

In the next section, the observations and data reduc¬ 
tion are described. Section 3 outlines the modeling of the 
observed spectra, and the resulting properties are discussed 
in Section 4. These are used to determine the remaining life¬ 
time in Section 5. Finally, we conclude on the nature of the 
WO stars in Section 6. 


2. Observations and data reduction 

All observations presented in this work have been ob¬ 
tained at the European Southern Observatory using the X- 
Shooter instrument (Vernet et al. 2011) on the Very Large 
Telescope. WR142 and WR30a were observed under pro¬ 
gram ID 093.D-0591, and all the other stars as part of 
the NOVA program for guaranteed time observations un¬ 
der program ID 091.C-0934. X-Shooter covers a wavelength 
range from 3 000 A to 25 000 A by directing the light in 
three separate arms: the UVB (3 000-5 500 A), VIS (5 500- 
10 000 A), and NIR (10 000-25 000 A) arms. 

To prevent detector saturation by the strong emission 
lines, the observations were split up into several shorter ex¬ 
posures. An overview of the exposure times and slit widths 
used is given in Table 2. All observations were done in nod¬ 
ding mode with a nod throw of 5". The slit was oriented 
along the parallactic angle for all observations. 

The data of all stars have been reduced using the X- 
Shooter pipeline v2.2.0, which produces flux-calibrated ID 
spectra. The flux calibration is performed using observa¬ 
tions of the spectro-photometric standard stars listed in 
Table 2 taken during the same night. Two of the stars 
(WR30a and LH41-1042) had another bright object in the 
slit, and the standard reduction was not sufficient. For these 
stars, we first obtained the 2D spectrum for each nodding 
position separately, without sky subtraction. We then sub¬ 
tracted the sky background using a clean part of the slit. 

Because the atmospheric dispersion corrector was un¬ 
available at the time of observations, the traces of the two 
stars mentioned above were not at a constant position on 
the slit as a function of wavelength in the UVB spectra. 
We therefore extracted the ID spectra by integrating the 
flux at each wavelength over a sufficiently large part of the 
slit. While this introduces additional noise, the quality of 
the data is high enough for this to be negligible. In the VIS 
arm the traces of both stars were at a constant position on 
the slit. Here, we folded the spectra in the wavelength di¬ 
rection, and fitted a Gaussian to the stellar signal. The ID 
spectra were then subtracted using the parts of the slit cor¬ 
responding to ±3.5cr covered by the average cross-dispersed 
point-spread function. Figure 1 shows an example, in which 
the second object in the slit is located around pixel 50. 
While the median flux of this object over the full wave¬ 
length range is very low, it has a noticeable effect on the 
extracted spectrum if not taken into account, in particu- 
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Table 2. Overview of observations. 


ID 

MJD 

At start exp. 

rri a 

-t exp 

(s) 

Slit width“ 
(") 

Spec. Standard 

Airmass 

WR102 

56490.069 

4 X 120 

0.8, 0.9, 0.9 

GD153 

1.0 

WR142 

56783.371 

8 X 200/260/260 

0.8, 0.9, 0.9 

LTT7987 

2.3 

WR93b 

56489.138 

4 X 600 

0.8, 0.9, 0.9 

EG274 

1.0 

BAT99-123 

56409.978 

10 X 180 

0.8, 0.9, 0.9 

LTT3218 

1.7 

LH41-1042 

56522.351 

10 X 180 

0.8, 0.9, 0.9 

FEIGE-110 

1.8 

WR30a 

56771.212 

8 X 200/260/260 

0.8, 0.9, 0.9 

LTT7987 

1.6 

SK188 

56522.388 

8 X 40 

0.5, 0.7, 0.4 

EG274 

1.5 


Notes. If a single value is given it is the same for all three X-Shooter arms; if three values are given it is for the UVB, VIS 
and NIR arms, respectively. 


Table 3. Derived reddening and total-to-selective extinction. 


ID 

E{B - V) 

Rv 

WR102 

1.26 

3.10 

WR142 

1.72 

2.85 

WR93b 

1.94 

3.25 

BAT99-123 

0.19 

3.10 

LH41-1042 

— 

— 


lar at the shorter wavelengths. For both stars, the flux of 
the contaminating object was negligible in the NIR, and 
the ID spectra were obtained using the nodding reduction, 
as this facilitates a better correction for the copious telluric 
lines. The resulting fluxes connect well with the VIS spectra 
which indicates that the contribution of the second object 
can indeed be ignored. 

The flux-calibrated spectra of the stars were extinction 
corrected using the CCM extinction laws (Cardelli et al. 
1989; O’Donnell 1994). A CMFGEN model (see Section 3.2) 
of the corresponding metallicity was used as a template for 
the slope of the spectrum (see Figure 2 for an example). 
The value of the total-to-selective extinction Ry was only 
adjusted if a proper dereddening could not be achieved us¬ 
ing the average value of Ry = 3.1. The derived values for 
E{B — V) and Ry are listed in Table 3. The estimated 
uncertainty is 10% for Ry and 0.1 dex for E[B — V). 

The spectrum of one of the single WO stars, LH41-1042, 
shows a steeper slope, but does not display spectral features 
from a companion. Although a faint star is detected in a 
UV image at a small projected distance from LH41-1042 
(see Figure A.l), the spectrum could not be corrected by 
assuming a Rayleigh-Jeans contribution from this object. 
Instead, we have artificially corrected the flux to match the 
model WO slope for LMC metallicity. The drawback of this 
approach is that the reddening cannot be determined for 
this object, implying a larger uncertainty in the derived lu¬ 
minosity. The flux correction for LH41-1042 is described 
in detail in Appendix A. The extinction-corrected, flux- 
calibrated spectra of all WO stars are presented in Figures 3 
and 4. 

2.1. Spectral classification 

Two quantitative classification schemes exist for the WO 
subclasses. The first divides the WO class in five subclasses, 
ranging from WOl to W05 (Kingsburgh et al. 1995). This 
classification scheme is based on the equivalent width (EW) 



Fig. 1. Example of the extraction of the ID spectrum of WR30a 
for one of the nodding positions. The cross-dispersed profile is 
shown in black, the Gaussian fit in red. The blue dashed lines 
indicate the region that is extracted (corresponding to ±3.5 (t, 
with a the standard deviation of the fitted Gaussian profile). 

ratio of Ovi A3811-34 A to Civ A5801-12 A and Ovi 
A3811-34 A to Ov A5590 A. 

Crowther et al. (1998) introduced a classification 
scheme for WC and WO stars in which the WO class is 
divided into four subtypes, from WOl to W04, and con¬ 
nects to the WC class at WC4. This classification is based 
on the same EW ratios as the scheme from Kingsburgh 
et al. (1995), but also includes the full width at half maxi¬ 
mum (FWHM) of Civ A5801-12 and the EW ratio of Ovii 
A5670 to Ov A5590. 

We adopt the classification scheme of Crowther et al. 
(1998) in this paper. The spectral types that we derive are 
given in Table 1 (see also Appendix B). 

3. Spectroscopic analysis 

First, we briefly discuss the morphological properties of the 
flux-calibrated spectra. Following this, we perform a de¬ 
tailed quantitative spectroscopic analysis. 

3.1. Morphological properties 

The WO spectra in Figures 3 and 4 show clear trends with 
spectral type and metallicity. These trends reflect changes 
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Fig. 3. Dereddened, flux-calibrated X-Shooter spectra in the UVB and VIS range (3 000— 10 000 A). The flux has been multiplied 
by an arbitrary factor for plotting purposes. Nebular emission and residuals from the sky subtraction have been clipped, and 
wavelength ranges affected by strong atmospheric features are displayed within the shaded areas. The spectrum of DRl has been 
rebinned to 1 A. 



log{Wavelength) (A) 


Fig. 4. As Figure 3, but for the NIR range (10 000 — 25 000 A). The spectra of BAT99-123, LH41-1042 and Skl88 have been binned 
to 1 A. The spectrum of DRl has been binned to 2 A, and only extends to 20 000 A due to the use of the X-Shooter A-band 
blocking filter. 
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log(Wavelength) (A) 

Fig. 2. Example of the extinction correction procedure for the 
highly reddened WR93b. The flux-calibrated spectrum before 
and after the extinction correction is displayed in red and black, 
respectively. The scaled continuum of a model spectrum is indi¬ 
cated by the dashed black line. 
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Fig. 5. Comparison of the primary and secondary spectral clas¬ 
sification criteria (x and y axis, respectively). The shaded areas 
indicate the regions where the two criteria agree on the spectral 
type. If these are not in agreement, the classihcation depends on 
the FWHM of Civ A5801-12 A. This is only the case for DRl, 
which we assign the spectral type W03, in agreement with pre¬ 
vious classifications. The dashed shaded region indicates where 
the WC sequence begins. 

in the physical properties of the star and in the region in the 
wind where the lines are formed, and provide information 
for the subsequent modeling (Section 3.2). 

Most notable is the increase in line width at higher 
metallicity. This is visible to a varying degree in all the 
spectral lines, and can be seen particularly well in the C iv 
and Hell lines in the NIR. The broader lines at higher 
metallicity reflect an increase in the terminal velocity of 
the outflows from the stars. 

The morphology of the characterizing Ovi A38I1-34 A 
emission also changes with metallicity. At low metallicities 
(IC 1613 and LMC, but also for the binary in the SMC) this 
doublet is clearly double-peaked, while it is fully blended 
in the spectra of the galactic stars. The O v A5590 A line 



Fig. 6. Behavior of the Hell, Ov, and Civ blend at 6560 A for 
different values of the temperature and the oxygen abundance. 
The observed spectrum is from WR142. 


profiles change from a roughly parabolic shape in the stars 
with a sub-galactic metallicity to a broad flat-topped shape 
in the MW stars. This reflects a change in the optical depth 
of the line-forming region: the line is formed in the optically 
thick region of the outflow for the lower metallicities, while 
for the galactic stars it is formed in optically thin regions. 

3.2. Modeling 

To perform a homogeneous quantitative spectroscopic anal¬ 
ysis of the WO stars we employ the CMFGEN code of Hillier 
& Miller (1998). This code iteratively solves the transfer 
equation in the co-moving frame, and accounts for effects 
such as clumping and line blanketing. Our fitting strategy 
has been described in detail in Paper I. Here we only report 
on relevant assumptions and changes in the applied diag¬ 
nostics. Our modeling approach aims to reproduce all the 
observed trends in the sample, as well as to provide a good 
fit to each of the individual spectra. 

In Paper I, the weak optically thin Hell line at 4859 
A was used as a diagnostic for the stellar temperature. 
This line is not recognizably present in the spectra of the 
LMC and MW stars that are analyzed here. Instead, we 
use the shape of the Heii A6560 A as a temperature probe. 
The blue wing of this line has contributions from O v and 
C IV, and the shape of the line profile can only be fitted by 
models with the correct combination of temperature and 
carbon and oxygen abundances. Together with the other 
abundance diagnostics, this allows us to constrain the tem¬ 
perature with an accuracy of about 20 kK (see Figure 6). 

To determine the luminosities of the stars, we adopt a 
distance of 50.12 kpc to the LMC (Gibson 2000), and dis¬ 
tances of 4.6, 1.75, and 3.4 kpc for WO-MW-1, 2, and 3, 
respectively (Drew et al. 2004). We derive spectrophoto- 
metric V-band magnitudes for the observed spectrum and 
compare these to the Simbad values. This gives an esti¬ 
mate of the flux loss and error in the flux calibration of the 
observations. We model the luminosity by matching the 
distance-corrected model flux to the magnitude-corrected, 
dereddened observed flux in the V-band wavelength region. 
The model flux is not only determined by the luminosity. 


5 




























F. Tramper et al.: Massive stars on the verge of exploding: the properties of oxygen sequence Wolf-Rayet stars 


Table 4. Properties of the single WO stars. 


ID 

Typical uncertainty 

logL 

[Lq) 

T* 

(kK) 

20 

R* 

{Rq) 

No 

Nhs 

10% 

No 

NHe 

10% 

Voo 

(km s“^) 
200 

logM 
(M0 yr-l) 
0.1 

log Rt 
(Rq) 

V 

logQo 

(s-l) 

logQi 

(s-l) 

logQ2 

(s-l) 

WR102 

WR142 

WR93b 

BAT99-123 

LH41-1042 

DRl 

K 4K+0.15 
o.^O_Q 23 

r oq+ 0.15 
^"^^-0.23 
c; Qn+0-15 

O.OU_Q 23 

5.20^°“® 

5-26l°o:l4 

5.68l°:S®r 

210 

200 

160 

170 

150 

150 

0.39 

0.40 

0.58 

0.47 

0.62 

1.06 

1.50 

1.00 

0.60 

0.63 

0.90 

0.35 

0.45 

0.16 

0.15 

0.13 

0.20 

0.06 

5000 

4900 

5000 

3300 

3500 

2750 

-4.92 

-4.94 

-5.00 

-5.14 

-5.05 

-4.76 

0.67 

0.79 

0.97 

0.82 

0.83 

0.46 

10.3 

11.3 

12.3 

7.3 
8.6 

5.0 

49.0 

49.0 

48.9 

48.9 

48.9 

49.5 

48.9 

48.8 

48.8 

48.7 

48.7 

49.3 

48.4 

48.3 

47.3 

47.6 

45.5 

48.0 


but is also affected by the mass-loss rate and stellar temper¬ 
ature. Therefore, like the determination of all other param¬ 
eters, the luminosity determination is an iterative process. 

Table 4 presents an overview of the best-fit parameters 
for all the single WO stars in our sample. A formal error 
determination is not possible, as we cannot perform a full 
exploration of parameter space. Table 4 gives an estima¬ 
tion of the uncertainty in the derived values, which reflect 
the sensitivity of the model spectra to changes in that pa¬ 
rameter. To estimate the uncertainty in the luminosity, we 
assume an error of 20% and 5% in the distance for the 
Galactic and extra-Galactic stars, respectively. Based on 
the comparison with the V-band magnitudes, we estimate 
an uncertainty of 10% in the flux calibration. The resulting 
error bars are given in Table 4. The error on the luminos¬ 
ity of LH41-1042 has been multiplied by two to reflect the 
absence of the extinction correction (see above). 

Figures showing the best-fit model for each star are 
given in Appendix C. Table 4 also gives values for the trans¬ 
formed radius i?*, the wind efhciency rj, and the ionizing 
fluxes (5 o,i, 2- The ionizing fluxes indicate the number of 
photons per second that are available to ionize hydrogen 
(Qo)) and singly (Qi) and doubly {Q 2 ) ionize helium. 

The transformed radius, introduced by Schmutz et al. 
(1989), is defined as: 


Rt — R* 


'^OD 

2500 kms"^ 


M 


2/s 


V/clO-^Moyr-l 


( 1 ) 


where the volume filling factor /c was first incorporated by 
Hamann & Koesterke (1998). The temperature T* is the 
temperature at radius i?* at the base of the wind, where 
the underlying regions can be assumed to be in hydrostatic 
equilibrium. Combinations of parameters that keep Rt con¬ 
stant produce very similar spectra, making it a very useful 
quantity. For instance, it implies that when Voo and i?* are 
well constrained, the line flux is determined by the value of 
M/^/Jc- In Section 5 we use this property to constrain /c. 

The wind efficiency parameter gives the ratio of the 
wind momentum once the flow has reached the terminal 
velocity (Mvoo) and the photon momentum (Ljc): 

, = ( 2 ) 

Thus, rj indicates the average amount of scatterings that 
photons undergo in order to drive the wind. The values of 
rj for the WO stars are in the range of 5 ~ 12, very similar 
to values found for WC stars (Sander et al. 2012). Multiple 
photon scatterings are expected for optically thick winds 


(de Koter et al. 1997; Vink & Grafener 2012; Grafener & 
Vink 2013). 

Overall, the observed spectra are well reproduced by 
our models, and allow us to constrain the temperature, 
surface abundances, and wind properties of the WO stars. 
Nevertheless, some of the observed spectral features are not 
fully reproduced. The Ovi A3811-34 A cannot be repro¬ 
duced while simultaneously htting the overall spectrum, 
and the flux in this line is underpredicted by a factor of 
~ 3 in our models. The cause of this is likely to be found 
in the susceptibility of the population of the upper level of 
the transition to X-ray excitation (see Paper I). Still, the 
strength of this line in our models does follow the observed 
trend with spectral type. 

The observed spectra show some emission lines that are 
not in the model spectra. These lines belong to transitions 
of the higher ionization stages of oxygen and carbon (i.e., 
Ovii, Oviii, and Gv), for which we do not include atomic 
models. These ionization stages are treated as auxiliary lev¬ 
els, to ensure that the populations of the levels that are 
modeled are accurate. While it is possible to include these 
high ionization stages in our models, their level populations 
are highly dependent on the soft X-ray radiation field, and 
thus the presence of, e.g., shocks. This should not affect the 
derived properties of the stars, but does result in a slightly 
higher uncertainty in the temperature. 

We found that for models with very high temperature 
(T* ^ 190 kK) the ratios between the Oiv A3404-12 A, 
O V A5598 A, and O vi A5290 A lines can no longer be re¬ 
produced. For the stars that have such high temperatures 
we determine the oxygen abundance by adopting models 
where the oxygen line ratios are closest to the observed 
values. The uncertainty in the obtained abundance is in¬ 
vestigated by fitting each of the individual oxygen lines, 
while adjusting the other parameters to preserve the over¬ 
all fit to the spectrum. The range in oxygen abundances 
that results from this approach is small (< 5% changes in 
the derived mass fractions), and thus we expect them to be 
accurate. 

Lastly, all spectra show to some extent a very broad 
emission feature at ~ 21000 A. This feature is a blend of 
several emission lines of G iv, G ill and He ii, yet almost 
none of our models show emission in this region, with the 
model for LH41-1042 being the exception. 


4. The properties of the single WOs 

In this section we discuss the derived properties of the WO 
stars and we compare our results to previous research. We 
then place the WO stars in the Hertzsprung-Russell dia- 



F. Tramper et al.: Massive stars on the verge of exploding: the properties of oxygen sequence Wolf-Rayet stars 


Table 5. Comparison of estimates of the terminal wind velocity. Table 6. Comparison of estimates of the surface abundances. 


ID 

'1^00 

(km s“^) 

Reference 

ID 

Nc/Nue 

No/Nue 

Reference 

WR102 

5000 

This work 

WR102 

1.50 

0.45 

This work 


5000 

Sander et al. (2012) 


0.51 

0.11 

Kingsburgh et al. (1995) 


4600 

Kingsburgh et al. (1995) 


0.44“ 

0.25“ 

Sander et al. (2012) 

WR142 

4900 

This work 

WR142 

1.00 

0.16 

This work 


5000 

Sander et al. (2012) 


0.52 

0.10 

Kingsburgh et al. (1995) 


5500 

Kingsburgh et al. (1995) 


0.44“ 

0.25“ 

Sander et al. (2012) 

WR93b 

5000 

This work 

WR93b 

0.60 

0.15 

This work 


5750 

Drew et al. (2004) 


0.95 

0.13 

Drew et al. (2004) 

BAT99-123 

3300 

This work 

BAT99-123 

0.63 

0.13 

This work 


4100 

Crowther et al. (2002) 


0.51 

0.11 

Kingsburgh et al. (1995) 


4300 

Kingsburgh et al. (1995) 


0.70 

0.15 

Crowther et al. (2002) 

DRl 

2750 

Tramper et al. (2013) 

DRl 

0.35 

0.06 

Tramper et al. (2013) 


2850 

Kingsburgh & Barlow (1995) 


0.63 

0.27 

Kingsburgh & Barlow (1995) 


Notes. Fixed value (not fitted). 


gram (HRD) and compare their position with those of the 
WC stars and with evolutionary tracks. Lastly, we discuss 
the mass-loss properties of the sample. 

4.1. Comparison with previous results 

Stellar temperatures have been derived for WR102 and 
WR142 by Sander et al. (2012), who found 200 kK for both 
stars, in agreement with our results. Crowther et al. (2000) 
report T* = 150 kK for BAT99-123 based on their model¬ 
ing of the far-UV to visible spectrum. The latter tempera¬ 
ture is based on a compromise between the UV and optical 
diagnostics, as the authors cannot simultaneously fit the 
O VI 1032-38 A resonance line (requiring 120 kK) and C iv 
7700 A (requiring 170 kK, in agreement with our value). 
To achieve a fit of both lines likely requires a detailed mod¬ 
eling of the wind acceleration and the clumping properties 
throughout the outflow. 

For all stars except LH41-1042 one or more measure¬ 
ments of the terminal wind velocity are reported in the 
literature (see Table 5). In most cases our values are lower 
than those previously reported and that are based on full- 
width at zero line emission and blue edge absorption of ul¬ 
traviolet P-Cygni profiles. Disparities compared to the first 
method may be related to a calibration issue. Those with 
the second method may point to a difference in the treat¬ 
ment of the velocity stratification of the outflow. Several 
studies report evidence for the presence of two acceleration 
zones in early WC star winds, in line with theoretical con¬ 
siderations (Schmutz 1997). Such a stratification may be 
modeled using a double-/3 law (e.g., Crowther et al. 2000, 
2002; Grafener & Hamann 2005) and can explain the higher 
Vao values from P-Cygni profiles that originate further out 
in the wind than the optical and near-infrared recombina¬ 
tion lines used by us. 

Various values of the surface abundances of WO stars 
have been reported in the literature, an overview of which 
is presented in Table 6. Most of these values are derived 
based on a comparison of equivalent width measurements 
to recombination line theory, and deviate from our results 
considerably. This may simply reflect that measuring the 
equivalent widths of WR lines is difficult due to line blend¬ 
ing and the poorly defined continuum. However, the val¬ 
ues derived by detailed atmospheric modelling by Crowther 



Fig. 7. Location of the single WO stars in the Hertzsprung- 
Russell diagram. Also indicated are the WC stars analyzed by 
Sander et al. (2012) and evolutionary tracks for solar metallicity 
and an initial rotational veloctity of 40% critical from Ekstrom 
et al. (2012). The helium zero-age main sequence (He-ZAMS) 
from BEC models (which include envelope inflation) is indicated 
for solar metallicity (solid line) and SMC metallicity (dashed 
line). 

et al. (2002) do agree with our values within the estimated 
10% uncertainty. 

4.2. Hertzsprung-Russell diagram 

Helium-burning stars that show the products of helium 
burning in their spectra are expected to be located very 
close to the helium zero-age main sequence (He-ZAMS) for 
most of the helium-burning lifetime. They only evolve to¬ 
ward hotter regions after exhausting the helium in their 
core (see, e.g.. Figure 11). 

Figure 7 shows the HRD with our results for the WO 
stars and the results from Sander et al. (2012) for Galactic 
WG stars. All WO stars are located on the hot side of the 
He-ZAMS, which is a first indication that they may be post¬ 
helium burning objects. The WC stars are located between 
the ZAMS and He-ZAMS, and thus appear to be too cold 
for their core-helium burning state. In the past, this dis- 
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Fig. 8. Mass loss versus luminosity relation with the location of 
the WO stars. Results for Galactic WC stars from Sander et al. 
(2012) and LMC WC stars from Crowther et al. (2002) are also 
indicated. Also plotted are the prediction for WC abundances 
of Nugis & Lamers (2000) with an additional metallicity 

dependence (Vink & de Koter 2005). The grey symbols indicate 
the Nugis & Lamers (2000) predictions for the WO stars, which 
underpredict the observed mass-loss rates. 


crepancy has often been attributed to the extended pho¬ 
tospheres of WR stars. However, the problem remains in 
recent studies that fully take this effect into account (e.g., 
Sander et al. 2012; Hainich et al. 2014). 

A possible solution to explain the positions of WC 
stars in the HRD is stellar envelope inflation near the 
Eddington limit. In this context, inflation refers to the ex¬ 
tended low-density envelopes in stellar models that reach 
the Eddington luminosity in their outer layers. This effect is 
predicted to be strongest at high metallicities (compare the 
solar and SMC metallicity He-ZAMS in Eigure 7). Grafener 
et al. (2012) could bring the predicted WR radii in agree¬ 
ment with the observations if the inflated sub-surface layers 
are clumped. Grafener & Vink (2013) discussed the solution 
topology of inflated envelopes with optically thick winds 
and concluded that there are two types of possible solu¬ 
tions: cool stars with clumped inflated envelopes, and hot 
stars without envelope inflation but a lower mass-loss rate. 
While the former complies with the properties of late-type 
WG stars, the latter may represent the small groups of WO 
stars and early-type WG stars. 

The evolutionary tracks for solar metallicity from 
Ekstrdm et al. (2012) are also plotted in Figure 7. The 
tracks do not reach the region of the HRD where most 
of the WO stars are located, which indicates that single¬ 
star models cannot currently explain them. However, if 
the mass loss of stars with masses less than ^ 32Mq is 
higher than used in the tracks, these stars might reach the 
WO star region. In this case, the WO stars seem to be 
the descendants of stars with initial masses in the range 
25Mq ^ Mini ^ 60Mq. The location of the He-ZAMS (for 
and Z — 0.2Zq) that is plotted in Figure 7 cor¬ 
responds to the helium-burning models used in Section 5. 
Towards higher luminosities, the He-ZAMS bends towards 
lower temperatures which is caused by inflation. Again, the 
location of the WO stars favors them to be post-core helium 
burning, regardless of any uncertainties in the metallicity. 



Fig. 9. Same as Figure 8, but with the WRl relation of Yoon 
& Langer (2005) adjusted to have a metallicity dependence of 
(Vink & de Koter 2005). 


The evolutionary state of the WO stars is further dis¬ 
cussed in Sections 5 and 6. 

4.3. Mass-loss properties 

Figure 8 compares the mass-loss rates found for the WO 
stars to those of the WG stars from Sander et al. (2012, 
MW stars) and Growther et al. (2002, LMG stars). Also 
plotted are various forms of the empirical mass-loss predic¬ 
tions from Nugis & Lamers (2000) for hydrogen-free WR 
stars: 

logM = —11 -I-1.29 log 3- 1.71ogXHe + 0.5 log Z. (3) 
Lq 

Here, Xfje is the helium mass-fraction, and Z = \ — Aine- 
Z is thus almost equal to the sum of the carbon and 
oxygen mass fractions. Sander et al. (2012) find that the 
mass-loss rates for their WG stars are compatible with 
these predictions when using their average carbon and oxy¬ 
gen mass-fractions (Xhs = 0.55, Z = 0.45, upper grey 
line in Figure 8). Crowther et al. (2002) find an offset of 
approximately —0.2 dex for the mass-loss rates of LMC 
WC stars, corresponding to a metallicity dependence of 
~ Zpg®. Vink & de Koter (2005) predict a scaling of the 
mass-loss rates with Zp^^ for WC stars with metallicities 
0.1 ^ Zpe/Zpe,© ^ 1, which we implement in Figure 8. 

The Nugis & Lamers (2000) rates with the Vink & de 
Koter (2005) metallicity scaling match the observed WC 
mass-loss rates from Sander et al. (2012) and Crowther 
et al. (2002) well. However, they severely underpredict the 
mass loss of the WO stars, which have a very low helium 
abundance. This is most notable for WR102, for which the 
predicted mass-loss rate is about 0.5 dex lower than the 
observed rate. This needs to be taken into account when 
comparing to evolutionary predictions that use the Nugis 
& Lamers (2000) mass-loss prescription. 

Figure 9 shows the same comparison as Figure 8, but 
with the mass-loss relation WRl from Yoon & Langer 
(2005), which is a scaled down version of the results of 
Hamann et al. (1995). We adjust this relation to have the 
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Fig. 10. Scaling of the terminal wind velocity in WO stars with 
metallicity. 


same metallicity dependence as applied in Figure 8, yield¬ 
ing: 

T >7 

log M = -12.73 -h 1.5 log — -h 0.66 log —^. (4) 

Lq ^Fe,0 

Again, the mass-loss rates of the WC stars are well rep¬ 
resented by this relation. The LMC and IC1613 WO stars 
also have a mass loss that is close to the predicted value. 
Only the mass-loss rate of the galactic WO stars is not well 
reproduced by Equation 4. 

In line with earlier studies (e.g., Crowther & Hadfield 
2006), the terminal wind velocities that we derive scale with 
metalliciy. This behavior is shown in Figure 10. The metal¬ 
licity dependence can be well represented by a linear rela¬ 
tion. 


5. Evolutionary state: remaining lifetime and final 

fate 

The derived properties of the WO stars allow us to estimate 
the remaining lifetime of these stars, as well as to predict 
the type of SNe they will produce. As helium burns at a 
more or less constant rate, the currently observed surface 
composition corresponds to material that was in the (fully 
mixed) convective core of the star at a time that is roughly 
proportional to (1 — Aue) x rue, with rue the duration of 
core-He burning. In Paper 1, we used this relation to put a 
limit on the current stage of helium-burning for DRl. 

In this work, we take a more detailed approach. We use 
the Binary Evolution Code (BEC, Yoon et al. 2006, 2010; 
Brott et al. 2011) to model the WR stars as non-rotating 
hydrogen-free helium stars. BEC is a state-of-the-art, one¬ 
dimensional hydrodynamic implicit Lagrangian code, and 
is well-suited to investigate stars that evolve close to the 
Eddingon limit (Kohler et al. 2015). 

In particular, convection is treated within the frame¬ 
work of the standard Mixing Length Theory (Bohm-Vitense 
1958) with the mixing length set to 1.5 times the local pres¬ 
sure scale-height. All the models were set to hydrostatic 
equilibrium. The models are evolved from the He-ZAMS 
until at least core-oxygen ignition. 



Fig. 11. Hertzsprung-Russell diagram with evolutionary models 
for WR102, which start with a helium star. The best-fit param¬ 
eters for WR102 are indicated. The different tracks correspond 
to different mass-loss rates. The dots indicate 1000 year time 
steps. 

We calculate evolutionary models for each of the WO 
stars, with the aim to reproduce the stellar temperature and 
luminosity at the point where the surface helium abundance 
of the models equals the observed values (see Table 7). We 
use the metallicities listed in Table 1. As shown before, 
the mass-loss rates from Nugis & Tamers (2000) that are 
normally used in these evolutionary models are not repre¬ 
sentative for our stars. We mitigate this by adopting a con¬ 
stant mass-loss rate throughout the evolution. We explore 
combinations of the mass-loss rate and initial helium star 
mass to reproduce the observed surface helium abundance, 
temperature and luminosity. 

Figure 11 shows the stellar evolutionary tracks for 
WR102 in the HRD, computed with an initial helium mass 
of Mne.ini = 22Mq and three values for the mass-loss rate. 
These were the best-fit models for the observed helium mass 
fraction, luminosity and temperature. Other combinations 
of initial helium-star masses and mass-loss rates were ex¬ 
plored but they were not in as good agreement with the 
observed stellar parameters. The surface mass fractions of 
helium (normalized to Ape + Ac -I-Aq = 1) during the late 
evolutionary stages are indicated along the tracks. Models 
for the other WO stars are shown in Appendix D. 

The models spend most of their lifetime close to the 
He-ZAMS. They lose mass as a consequence of the ap¬ 
plied mass-loss rate causing the luminosity to drop during 
the core-helium burning phase. After core-helium exhaus¬ 
tion, the stars contract on their thermal timescale and be¬ 
come hotter and brighter. The models are hardly inflated 
at this stage of evolution because the iron opacity peak (at 
T ^ 200 kK) which is responsible for inflation is only par¬ 
tially present in the stars. At some point during the overall 
contraction phase the helium-shell ignites and the tracks 
eventually turn redwards because of the mirror principle 
(Kippenhahn & Weigert 1990). When the models become 
cooler, envelope inflation may again play a role. The mod¬ 
els were computed until oxygen burning and the star will 
not change its position in the HRD significantly after the 
models stop, as the envelope structure does not change any 
further until the star explodes as a supernova. 
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Table 7. Surface mass fractions, estimates of helium star mass and remaining lifetime, and constraints on the mass loss parameters. 


ID 

Vhb 

Ac 

Ao 

AfHe.ini 

(Me) 

Afflnal 

(Mq) 

tSN 

(yr) 

Mevol 

(Moyr-l) 

/c 

WR102 

0.14 

0.62 

0.24 

22.0 

9.8 

1500 

2.8 X 10“® 

> 0.4 

WR142 

0.26 

0.54 

0.21 

17.0 

8.8 

2 000 

1.7 X 10“® 

> 0.4 

WR93b 

0.29 

0.53 

0.18 

17.0 

8.8 

8 000 

1.7 X 10“® 

> 0.3 

BAT99-123 

0.30 

0.55 

0.15 

15.0 

7.7 

7000 

1.4 X 10“® 

> 0.3 

LH41-1042 

0.22 

0.60 

0.18 

17.0 

8.4 

9 000 

1.8 X 10"® 

> 0.4 

DRl 

0.44 

0.46 

0.10 

23.0 

15.4 

17000 

1.8 X 10"® 

0.1 


The stellar tracks in Figure 11 that reproduce the ob¬ 
served stellar parameters do so after the models have ex¬ 
hausted the helium in their core (see also hanger et al. 
1988; hanger 1989). In other words, according to our mod¬ 
els WR102 is a post-core helium burning star and has a 
remaining lifetime of less than 2000 years. Furthermore, 
tracks with mass-loss rates lower than ~ 2.7 x lO“®M 0 yr“^ 
never reveal the layer that corresponds to the observed he¬ 
lium mass-fraction. Combined with the determined value of 
M / vTc, this constrains the volume filling factor to /c > 0.4. 

As the migration to higher temperatures after core¬ 
helium exhaustion occurs on a very short timescale, the 
surface abundances do not change significantly during this 
period. This implies that a few WC stars with similar sur¬ 
face abundances may be expected to exist. Koesterke & 
Hamann (1995) indeed find carbon surface mass fractions 
for WC stars in the range 0.2 —0.6, i.e. covering the 0.4 —0.6 
range that we find for the WO stars (Table 7). The oxy¬ 
gen abundance is not well determined for WC stars (e.g., 
Crowther 2007) so a direct comparison with our results can¬ 
not be made. 

Table 7 gives the initial helium star mass, predicted 
remaining lifetime, and constraints on the clumping fac¬ 
tor based on the calculated evolutionary models for all the 
stars. According to our models, all WO stars except maybe 
DRl (which is close to the helium terminal-age main se¬ 
quence) are in the core-contraction phase after core-helium 
burning. They are likely to explode in less than 10"^ years. 
The initial helium star masses are in the range Mue.ini ~ 
15—25 Mq. This corresponds to the helium-core masses pre¬ 
dicted for stars with an initial mass of Mi„i ~ 40 — 60 Mq 
(Ekstrom et al. 2012). However, only the models with very 
rapid rotation (40% of critical. Figure 7) reach the He- 
ZAMS and hotter regions, while non-rotating models are 
still covered by a hydrogen-rich envelope. As most massive 
stars are only modestly rotating (e.g., Ramfrez-Agudelo 
et al. 2013), this may point to a much higher mass loss 
(either through the stellar wind or binary interactions) in 
stages prior to the WO phase. 

Groh et al. (2014) modelled the spectrum of a non¬ 
rotating star with an initial mass of 60 Mq during various 
stages of its evolution, based on the parameters predicted 
by the Ekstrom et al. (2012) evolutionary tracks. These 
tracks produce a WO spectrum during a very short stage 
prior to core-helium exhaustion. After core-helium exhaus¬ 
tion the tracks again enter a longer WO phase that lasts 
until the star explodes. We find that the WO star with an 
estimated initial mass close to 60 Mq (WR102) has a lu¬ 
minosity that is much lower then those predicted by the 
Ekstrom et al. (2012) tracks (see Figure 7). This indicates 
that more mass has been lost during its evolution compared 
to the theoretical predictions. 


While all observed WO stars are post-core helium burn¬ 
ing, but this does not exclude that some of the WC stars 
are in this phase as well, and thus could explode without 
ever becoming a WO star. Whether the whole star contracts 
after core-helium exhaustion (i.e. becomes a WO star) is de¬ 
pendent on the mass loss prior to, and during, the helium¬ 
burning phase and must be much higher then the currently 
used Nugis & Earners (2000) predictions. 

When they end their lives, most of the WO stars have 
predicted masses below \QMq. The supernovae will almost 
certainly be of type Ic, as the current fractions of helium 
at the surface are already too low to produce type Ib SNe 
(Dessart et al. 2011). Although Sander et al. (2012) report 
on a very high rotational velocity for WR102 and WR142, 
we do not find any indication of rapid rotation for the WO 
stars. We can therefore not conclude on the possibility of 
the production of long-duration GRBs during core-collapse. 

5.1. Progenitors of 14/0 stars 

As shown above, the likely progenitors of WO stars have 
an initial mass of Mini = 40 — 60Mq . This is in very good 
agreement with the current view of massive star evolu¬ 
tion, for instance as summarized by Langer (2012). The 
evolutionary sequences proposed by Gonti (1975) and later 
refined by Langer et al. (1994) and Sander et al. (2012) 
also predict WO stars to be the descendants of stars with 
Mini = 45 - 6 OM 0 . 

As mentioned before, single-star evolutionary tracks 
that start from the ZAMS cannot reproduce most of the 
WO stars. However, the evolutionary path of massive stars 
and the mass-loss rates at each of the various stages of evo¬ 
lution are highly uncertain. In particular, whether a star 
undergoes an LBV phase, and how much mass is lost during 
that phase, is currently poorly understood. The question is 
therefore how the naked helium stars that are needed to 
explain the WO phenomena are formed. While some WO 
stars may indeed originate from single-star evolution, bi¬ 
nary interactions may also play a role (Sana et al. 2012). 

If the WO star progenitor is initially the most mas¬ 
sive star in a binary system, it might strip its envelope by 
transferring mass to a massive companion (e.g., Paczyhski 
1971). The hydrogen-rich material would rejuvenate the 
companion star, making it appear younger (e.g., Hellings 
1983; Dray & Tout 2007; de Mink et al. 2014). This sce¬ 
nario might be appropriate for the two WO stars that are 
in a short-period binary with an 0-star companion, WR30a 
and Skl 88 . The spectral signature of their companions are 
clearly visible in the observed spectra (see Figure 3). A 
quantitative analysis of these stars in a later work may re¬ 
veal the potential role of binarity in the evolution of WO 
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Fig. 12. Evolution of the surface mass fractions in the helium 
star model of WR102 since the onset of core-helium burning. 
The observed mass fractions of helium, carbon, and oxygen are 
indicated. The model over predicts the Xq/Xq ratio. 


star progenitors. The single WO stars that are discussed in 
this paper show no spectral signature of a massive compan¬ 
ion. However, we cannot exclude the presence of a fainter 
intermediate-mass companion as in the WO stars might 
have lost their hydrogen envelope in a common envelope 
interaction. In that case the companion star is expected to 
still orbit the WO star but would be too faint to be seen 
in the spectrum. If the WO stars do instead result from 
single-star evolution, an LBV phase after core-hydrogen ex¬ 
haustion may be the likely explanation to reveal the core at 
the onset of helium burning (as is the case in Groh et al. 
2014). 

5.2. The ^^C(q;, 7)^®0 thermonuclear cross-section 

Apart from constraining the lifetime of the WO stars, the 
observed surface abundances can be used to constrain the 
elusive ^^C(q;,7)^®0 thermonuclear reaction rate (Grafener 
et al. 1998). The currently used value (0.632 times that of 
Caughlan et al. 1985) is based on the solar abundance pat¬ 
tern between oxygen and calcium, but has an uncertainty 
of about 30% (Weaver & Woosley 1993). Deviations of this 
rate strongly influence the supernova yields as well as the 
pre-supernova evolution (Tur et al. 2007). With their sur¬ 
face abundances corresponding to the core abundances far 
into the helium-burning stage, WO stars offer a unique op¬ 
portunity to provide direct constraints on the nuclear re¬ 
action rate. Figure 12 shows the helium, carbon, and oxy¬ 
gen surface mass fractions as a function of time for the 
helium star model of WR102. When compared to the ob¬ 
served mass fractions, it is clear that the model overpre¬ 
dicts the oxygen-to-carbon ratio. The same effect is seen in 
the models for all the other stars. This indicates that the 
actual ^^C(q;, 7)^®0 reaction rate must be lower than the 
value that is currently being used. Placing more firm con¬ 
straints on this reaction rate will be the focus of a separate 
study. 


6. Summary and conclusions 

In this paper we have presented a detailed spectroscopic 
analysis of single WO stars. We have constrained the evo¬ 
lutionary status of the sources using tailored evolutionary 
models. For the stellar properties we find that: 

— WO stars are extremely hot, with temperatures ranging 
from 150 kK up to 210 kK. In the Hertzsprung-Russell 
diagram they are located at the hot side of the helium 
zero-age main sequence; 

— The luminosities range from 1.2x 10^ Lq to 3.2x 10® Lq, 
comparable to the luminosities of WC stars. Only DRl 
is much brighter than WC stars, with a luminosity of 
5 X 10® Tq. 

— The helium surface mass fraction is typically 20-30%, 
but ranges from 44% for the coolest star to as low as 14% 
for the hottest star. The oxygen mass fractions reach 
values up to 25%; 

— The stars lose mass at a rate of 0.8—1.5 x 10“®Mq yr“^, 
which is higher by a factor of about 2-3 than suggested 
by the predictions by Nugis & Tamers (2000); 

— The surface abundances of WO stars correspond to ma¬ 
terial that was formed very late into the core-helium 
burning stage, and the stars must have experienced 
severe mass loss during this stage to reveal the ob¬ 
served layers. This translates to strict constraints on 
the volume filling factor of the wind medium, requiring 
/c > 0.3 for most of the stars. 

The derived properties of the WO stars can be repro¬ 
duced for helium stars with Mne.ini 15 —25M0, and sug¬ 
gest initial masses of Mini ~ 40 — 60 Mq. Their extremely 
high temperatures are consistent with the contraction of the 
star after it exhausts its helium in the core. Together with 
the low surface helium abundances, this firmly establishes 
that WO stars are post-core helium burning massive stars. 
They will likely explode as type Ic supernovae in 10® — 10^ 
years. Thus, the WO class indeed represents a short, final 
stage in the evolution of massive stars. 
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Appendix A: Flux correction LH41-1042 



Fig. A.l. Archival HST/WFC3 image of LH41-1042 (F225W 
filter, proposal ID 12940, PI Massey). North is up and East to 
the left. The two X-Shooter slit positions are indicated. The 
inset shows a zoomed image of the boxed area around LH41- 
1042, and shows the nearby star that may be contaminating the 
spectrum. 

As discussed in Section 2, the slope of the spectrum of 
LH41-1042 before extinction correction is steeper than that 
of WO star models for an LMC metallicity. This may either 
be a result of a bad flux calibration, or a second object 
may be contaminating the spectrum. Inspection of a Hubble 
Space Telescope (HST) near-UV image shows that there is 
indeed a faint star very close to LH41-1042 (Figure A.l). 
We therefore first tried to correct the spectrum by assuming 
that this object is responsible for the steep slope. 

We assume that the flux contribution of the contam¬ 
inating object follows the Rayleigh-Jeans approximation 
{F oc A“^) at the X-Shooter wavelength range. This seems 
justified, as the slope of the uncorrected spectrum is steeper 
than the slope of the continuum of an LMC WO-model 
(which has a free-free emission component and is less-steep 
than a Rayleigh-Jeans slope) even before dereddening. 

As we cannot determine the reddening from our spec¬ 
trum, we use the average value for the region found by 
Massey et al. (2005) of Ay = 0.4 and a standard total-to- 
selective extinction of Ry = 3.1. We assume a luminosity of 
1.8 X 10® Lq for the WO star, equal to that of BAT99-123. 
While there is no a-priori reason for the two stars to have 
the same luminosity, the resulting model flux is in agree¬ 
ment with the dereddened flux in the near-IR region where 
the contribution of the contaminator becomes negligible. 

The correction is then done by testing different ratios 
of the flux contributed by the WO star and the contam¬ 
inator. Ideally, the flux ratio measured from photometry 
should be used. However, the only image where the indi¬ 
vidual stars are resolved is the HST/WFC3 image shown 
in Figure A.l, which uses a filter with an effective wave¬ 
length in the near-UV. From this image, the flux ratio is 
Fcont/Fwo = 0.1 at 2250A. As a much higher flux ratio is 
needed to recover the observed spectrum, this indicates that 



Fig. A.2. Determination of the flux contribution of the com¬ 
panion of LH41-1042. Plotted are the dereddened flux (black), 
the assumed LMC WO model continuum (blue), the derived 
Rayleigh-Jeans contribution (red), and the combined continuum 
from the WO and Rayleigh-Jeans contributions (green). 

the flux of the contaminating object peaks at a wavelength 
between 2250 A and the X-Shooter wavelengths (between 
~ 2500 — 3000 A), and no longer in the Rayleigh-Jeans tail. 
The contaminating object is therefore likely a B dwarf or 
giant. This means that the measured flux ratio can not be 
used to estimate the flux ratio in the X-Shooter wavelength 
range. Instead, we try out different combinations of flux ra¬ 
tios to obtain a combined spectrum of the WO model and 
a Rayleigh-Jeans contribution that matches the observed 
spectrum (Figure A.2.) 

Using this strategy, we can get a good representation of 
the slope of the observed spectrum (Figure A.2). However, 
when the Rayleigh-Jeans contribution is subtracted from 
the observed spectrum, the emission lines become roughly 
twice as strong as in any of the other WO stars. As oth¬ 
erwise the spectrum of LH41-I042 does not show unusual 
features, we conclude that this is unlikely to be physical. We 
therefore assume that the steep slope is a result of the flux 
calibration, and not of contamination by the faint nearby 
object. We correct for this by artificially altering the slope 
to the correct value. While this results in a much more re¬ 
alistic spectrum, the luminosity and mass-loss rate that is 
derived from the modeling are much more uncertain than 
those of the other WO stars. 

Appendix B: Equivalent widths and spectral 

classification 

Appendix C: Best-fit models 
Appendix D: Helium-burning models 
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Table B.l. Equivalent width {W\) measurements of the lines needed for the spectral classification. 


ID 

SpT 

Wa(Ovi A3811-34) 

(A) 

Wa(Ov A5590) 

(A) 

Wa(Civ A5801-12) 

(A) 

WR102 

W02 

1415 ±28 

196 ±4 

156 ± 1 

WR142 

W02 

1087±17 

199 ±2 

303 ±2 

WR93b 

W03 

475 ±5 

230 ±4 

1627 ± 7 

BAT99-123 

W03 

235 ±4 

104 ± 1 

2115 ±38 

LH41-1042 

W04 

171 ±4 

172 ±2 

1869 ± 28 

DRl 

W03 

146 ± 28 

125 ±5 

417 ± 11 



Fig. D.l. Same as Figure 11, but for WR142. 



Fig. D.3. Same as Figure 11, but for BAT99-123. 



Fig. D.2. Same as Figure 11, but for WR93b. 



Fig. D.4. Same as Figure 11, but for LH41-1042. 
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Fig. D.5. Same as Figure 11, but for DRl. 
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Fig. C.l. Best model of WR102. 




Fig. C.2. Best model of WR142. 
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Fig. C.3. Best model of WR93b. 




Fig. C.4. Best model of BAT99-123. 
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Fig. C.5. Best model of LH41-1042. 
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